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ABSTRACT 

We analyze the dynamical expansion of the H II region and outer photodissociation region (PDR) 
around a massive star by solving the UV and FUV radiation transfer and the thermal and chemical 
processes in a time-dependent hydrodynamics code. We focus on the physical structure of the shell 
swept up by the shock front (SF) preceding the ionization front (IF). After the IF reaches the initial 
Stromgren radius, the SF emerges in front of the IF and the geometrically thin shell bounded with the 
IF and the SF is formed. The gas density inside the shell is about 10 1-2 times as high as the ambient 
gas density. Initially the dissociation fronts expands faster than IF and the PDR is formed outside the 
H II region. Thereafter the IF and SF gradually overtakes the proceeding dissociation fronts (DFs), 
and eventually DFs are taken in the shell. The chemical composition within the shell is initially atomic, 
but hydrogen and carbon monoxide molecules are gradually formed. This is partly because the IF and 
SF overtake DFs and SF enters the molecular region, and partly because the reformation timescales of 
the molecules become shorter than the dynamical timescale. The gas shell becomes dominated by the 
molecular gas by the time of gravitational fragmentation, which agrees with some recent observations. 
A simple estimation of star formation rate in the shell can provide a significant star formation rate in 
our galaxy. 

Subject headings: Circumstellar matter - H II regions - ISM: molecules - STARS : formation 



1. INTRODUCTION 

As the H II region expands around a massive star, the 
shock front (SF) emerges preceding the ionization front 
(IF) . This SF sweeps up the ambient gas and the gas shell 
is formed at the edge of the H II region. Many authors 
have studied the scenario that the shell becomes unstable 
and next star formation is triggered around H II region 
(collect and collapse model, see, e.g., Elmegreen & Lada 
1977, Whitworth et al. 1994, Elmegreen 1998, and refer- 
ences therein). Recently, Deharveng et al.(2003, hereafter 
DV03) observed the fragmented molecular shell around the 
classical H II region, Shl04 and show the young stars (clus- 
ter) is formed in the core of one fragment. They argue that 
this is the evidence of collect and collapse model. Now we 
can refine the theory for the evolution of the shell and 
compare it with the detailed observation. Especially, we 
can focus on the amount of the fragmented molecular gas 
they observed. The SF initially emerges in front of the 
IF, where sufficient FUV photons to dissociate molecules 
are available. Therefore, it is not clear whether the SF 
can gather the molecular gas in the shell as observation 
indicates or not. 

Roger & Dewdney (1992) and Diaz-Miller, Franco & 
Shore (1998) studied the time-dependent expansion of H II 
and PDR solving the radiative transfer of UV and FUV 
photons. Although these works do not include hydrody- 
namics, they have shown that the IF gradually overtakes 
the dissociation front (DF) and this means that IF and 
preceding SF gradually enter the molecular region. The 
numerical study of the hydrodynamical expansion of the IF 
has been studied since 1960's (e.g. Mathews 1965, Lasker 
1966, Tenorio-Tagle 1976, Tenorio-Tagle 1979 , Franco et 
al. 1990 etc.). These works successfully show the various 
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dynamical aspects of the expansion in both homogeneous 
and inhomogeneous ambient medium. However, these do 
not include the outer photodissociation region (PDR) and 
thermal processes dominant in PDR. 

We perform the time-dependent calculation including 
the IF, DF, and the shell, which has been very limited. 
We solve the UV and FUV radiation transfer and hydro- 
dynamics numerically to investigate the structure and evo- 
lution of the shell and PDR as well as H II region. Since the 
compression rate behind the SF depends on the thermal 
processes, we include the thermal processes which domi- 
nate either in H II region or PDR. In this letter, we con- 
sider the simple situation that there is one massive star 
within the homogeneous ambient medium, which can be 
compared with Shl04 observed by DV03. The detailed 
quantitative aspects of the evolution will be given in subse- 
quent paper (Hosokawa & Inutsuka 2004, hereafter, Paper 
II). 

2. NUMERICAL METHODS 

We use a one-dimensional spherically-symmetrical nu- 
merical method. The numerical scheme for the hydrody- 
namics is based on the 2nd-order Lagrangian Godunov 
method (see, e.g., van Leer 1979). We use the on-the- 
spot approximation for UV and FUV radiation transfer. 
We assume that all hydrogen molecules are in the ground 
vibrational level of X 1 S g and ortho/para ratio is 3:1 (Diaz- 
Miller, Franco & Shore 1998). For the Lyman-Werner 
bands to photodissociate hydrogen molecule, we solve the 
frequency dependent transfer using the representative set 
of lines. The UV and FUV photon luminosity of the 
central star is adopted from Diaz-Miller, Franco & Shore 
(1998) for the case of Z = 1Z@. The main thermal pro- 
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cesses included in the energy equation is UV/FUV heat- 
ing (e.g., H photoionization, photoelectric heating) and ra- 
diative cooling (e.g., H recombination, Ly-a, 01 (63.0/im, 
63.1/xm) , Oil (37.29^m), CII (23.26^m, 157.7/zm), H 2 , 
and CO, see Hollenbach & McKee 1979, 1989, Koyama & 
Inutsuka 2000). Non-equilibrium reaction equations are 
implicitly solved for the species of e, H, H+, H 2 , C + , and 
CO. The ionization rate of O is assumed to be the same 
as that of H. We adopt the simple approximation for the 
dissociation process of CO molecule given by Nelson & 
Langer (1997). The hydrodynamic equations and the ra- 
diative transfer equation, energy equation, and reaction 
equations arc combined following Tenorio-Tagle (1976). 
The dust extinction is included only outside the IF and 
the dust temperature is calculated following Hollenbach, 
Takahashi & Ticlcns (1991). We have checked our numer- 
ical code with the well-studied problems (e.g., Franco et 
al. 1990). 

3. RESULTS OF THE NUMERICAL CALCULATION 

In this letter, we consider one typical case where there 
is one massive star of T cS = 40000 K (M — 41 M Q ) in 
the ambient medium. This corresponds to the 06V star, 
which is the central star of Shl04. The initial ambient 
number density of the hydrogen nucleon is unknown, and 
we adopt a typical value of the giant molecular cloud, 
tch.o = 10 3 cm~ 3 (riH 2 ,o = 500 cm -3 ). As in the stan- 
dard picture (e.g. Spitzer 1978), the IF expands rapidly 
as the weak R-type front in the recombination timcscale, 
ircc ~ 100 yr for njj.o = 10 3 cm~ 3 (formation phase). As 
the IF reaches Stromgren radius, R st = 0.56 pc, the pho- 
toionization rate and the recombination rate in the H II 
region becomes equal, the H II region begins to expand 
owing to the pressure difference between H II region and 
outer PDR or molecular cloud. At this phase, the SF ap- 
pears in front of the IF and the IF changes to the weak 
D-type front (expansion phase). The expansion law in this 
phase is given by 



Rw(t) = R gt 1 + 
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(e.g. Spitzer 1978), where Chii is the sound speed in the 
H II region. The dynamical timcscale is now given by 



yn 



Rst /Chii ~ 10 5 yr. Fig.l shows the time evolu- 
tion of the physical quantities of the gas in the expansion 
phase. The radius of Shl04 is 4 pc and the IF reaches 
this radius at t = 0.7 Myr in our calculation. For each 
snapshot, we can see the H II region (T ~ 10 4 K), the 
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the outer PDR (T ~ 100 K) and the outermost molecular 
cloud (T ~ 10 K). The SF emerges just in front of the 
IF. The density behind the SF significantly increases from 
the ambient value owing to the line cooling, and the gas 
density within the shell is about 10-100 times as high as 
the ambient density. The more detailed inner structure of 
the shell is shown below (Fig. 3). The geometrical thick- 
ness of the shell, h increases throughout the calculation, 
h ~ 0.02 - 0.2 pc. The PDR (T ~ 100 K) becomes narrow 
as the IF expands, that is, the IF gradually catches 




Fig. 1. — The snapshots of the gas-dynamical evolution. In each 
panel, five snapshots represents the profiles at i = 0.1,0.3,0.5,0.7, 
and 0.9 Myr respectively. 

up with the preceding DFs (e.g. Roger & Dewdney 1992, 
Diaz-Miller, Franco & Shore 1998). The upper panel of 
Fig. 2 represents the overtaking of the IF more clearly. The 
position of the SF is very close to that of the IF. As the 
system switches from the formation phase to the expansion 
phase, the timescale of the expansion of the IF suddenly 
becomes longer to idyn- Even in this phase, the preceding 
DFs continue to expand. Then the region between IF and 
DFs appears as a PDR. However, the DFs gradually slow 
down because the FUV flux at each DF decreases owing 
to the geometrical dilution and the dust extinction in the 
PDR. The IF expands because of the pressure difference 
between the H II region and the outer PDR, and the de- 
celeration is more slower according to eq.(l). 
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Fig. 2. — Upper panel : The time evolution of the various 
front positions. The solid and dotted lines mean the position of 
the IF and SF respectively. The broken (dot-solid) line represents 
the DF of H 2 (CO). To define the DFs of H 2 and CO, we use 
X H2 = 2n H2 /n Hnuc = 0.5 and X co = n co /n Crluc = 0.5. The 
shaded region corresponds to i < (Gp) -1 / 2 within the shell, and 
hence, the gravitational fragmentation is expected (see §4.1). Mid- 
dle panel : The time evolution of the column density of each region. 
Thin contour lines represent the position where Xjj 2 (broken), Xqq 
(dot-solid) = 0.2, 0.5 and 0.8. Lower panel : Same as the med- 
dle panel but for showing the gravitationally unstalbc region. The 
shaded region corresponds to the region where t < (Gp) -1 / 2 . 



The lower two panels of Fig. 2 represents the time evolution 
of the column density of each region. The column density 
of the H II region gradually decreases because the number 
density decreases as nun oc -R IF 3 ^ 2 in the expansion phase, 
where Rif is the radius of the H II region. The column 
density of the gas shell a increases. The ambient mass 
swept up by the time t is proportional to Rwit) 3 and the 
mass within the H II region is oc 7t.hii-Rif(£) 3 ^ RiF(t) 3 ^ 2 - 
Then the most of the swept-up mass does not flow into 
the H II region but remains in the shell. The calculated 
mass of the H II region is 430M Q and the shell mass is 
9OOOM0 at t = 0.7 Myr (including both PDR and molec- 
ular region), which is in good agreement with the obser- 
vation (see DV03). The PDR initially spreads outside the 
shell, but H 2 and CO molecules are gradually accumulated 
from the outer side of the shell. The accumulation of the 
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Fig. 3. — The distribution of gas density, pressure, temperature, 
and the chemical compositions at t = 0.3 Myr. In the lowest panel, 
Xh 2 and Xco are molecular ratios (see the caption of Fig. 2) and 
Xhii is the ionization rate defined as n H + /fiH nuc • 



SF. Therefore, the accumulation is not simply explained 
by the fact that the IF and the SF overtakes the DFs 
and that the SF enters the molecular region. The gas 
density within the shell is so high that the reformation 
timcscalcs of the molecules become short enough. The 
reformation timescale of the hydrogen molecule is about 
~ 10 6 yr for nn,o ~ 10 3 cm~ 3 . The dynamical timescale is 
about idyn ~ 10 5 yr, then the reformation of H 2 on the 
grain surface is not important with the ambient number 
density. However, the reformation timescale in the gas 
shell, <~ 10 4 ~ 5 yr, is shorter than td yn , therefore the for- 
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mation of H2 molecule is significant in the shell. 

Fig. 3 shows the detailed structure of the physical quan- 
tities and the chemical composition of the shell at t = 
0.3 Myr. As this figure shows, the molecular fraction 
in the shell is higher than the value in the PDR out- 
side the shell owing to the rapid reformation. As H II 
region expands, the DFs are taken in the shell and the 
denser (tih > 10 4 cm~ 3 ) PDR is formed in the inner side 
of the shell. Once the molecules are reformed in the 
shell, the FUV radiation is consumed by the photodissoci- 
ation of these molecules, which decelerates the expansion 
of the preceding DFs to accelerates the accumulation of 
the molecules. About 80 % (55 %) of the hydrogen (car- 
bon) atoms within the shell exist as H2 (CO) molecule at 
t = 0.7 Myr in our calculation. Fig. 3 also shows other in- 
teresting features of the shell. The thermal pressure inside 
the shell is high and about constant at the value of the H II 
region. The gas temperature has negative gradient from 
the inner side to the outer side of the shell, because the 
FUV radiation is attenuated owing to the dust extinction 
through the shell, and the photoelectric heating decreases. 
Conversely, the gas density has positive gradient to the 
outer side of the shell. 

4. DISCUSSION 

4.1. Dust in H II Region 

In this section, we discuss the possible effect of dust 
grains in the H II regions (e.g. Spitzer 1978). For ex- 
ample, the grains can be driven by the radiation pressure 
during the expansion of the H II region (e.g., Arthur et 
al. 2004). Arthur et al. have shown the small grains 
are destroyed by the sublimation but the large graphite 
grains can exist close to the star. In contrast, the recent 
observation by Roger et al. (2004) found no far-infrared 
component correlated with H II region, Shl70. The photo- 
electric heating works even in H II regions comparatively 
with the photoionization heating (Weingartner & Draine 
2001) and is sensitive to the uncertain abundance of the 
small dust grains. If we adopt the dust absorption cross 
section, 6 x 10~ 22 cm 2 (H nuclcon) -1 in H II region, which 
we adopt in the PDR and the outer region, the radius of 
the H II region becomes slightly smaller, but the differ- 
ence from the case without dust in H II region is always 
less than 5%. The temperature in the PDR decreases by 
10-20% because of the attenuation of FUV photons by 
dust grains in the H II region. The gas density within 
the shell rises by a factor, 1.4-1.5. Theses changes in the 
shell, however, only slightly promotes the accumulation 
of the molecular gas. If we include the standard photo- 
electric heating rate and dust recombination cooling rate 
in the H II region (Bakes & Tielens 1994), the tempera- 
ture near the central star in the H II region rises, but the 
temperature near the gas shell hardly changes. Therefore, 
the existence of dust grains do not significantly alter the 
results in this paper. 

4.2. Fragmentation of the Shell 

The triggered star formation scenario predicts that the 
shocked layer fragments and these fragments collapse to 
stars of the second generation. Although the fragmenta- 
tion of the shocked layer is studied by many authors, it is 
still uncertain which instability actually occurs and which 



induces the next star formation. Here, we assume that 
non-gravitational instabilities (e.g. decelerating shock in- 
stability (DSI)) themselves are not the essential ones for 
the sequential star formation. Elmegreen (1989) showed 
that the gravitational instability couples with DSI for the 
decelerating layer using the linear analysis, but MacLow & 
Norman (1993) pointed out that the DSI quickly saturates 
in the non-linear phase. 

However, Garcia-Segura & Franco (1996) showed the 
DSI is strongly modified with the presence of the IF. The 
shell rapidly fragments and the finger-like structures are 
formed in their 2D calculation. Furthermore, the shadow- 
ing instability and corrugation instability of the IF were 
studied by Williams (1999, 2002) and these instabilities 
also can disturb the IF and the shell. Although these 
complexities are beyond the scope of this paper, it will 
be interesting to study the effect of these instabilities on 
the evolution of the DF and the post-shock layer with 2D 
or 3D calculations in the future. 

Elmegreen & Elmegreen (1978), Miyama, Narita & 
Hayashi (1987), Lubow & Pringle (1993), and Nagai et 
al.(1998) have studied the stability of the dense gas layer 
and show that the earliest gravitational instability occurs 
with a wavelength which comparable to the layer thickness 
and the growth rate is ~ (Gp)^ 1 ^ 2 . The calculated shell 
suffers from this instability, and the shaded region in Fig. 2 
represents the unstable region, where t < (Gp)- 1 ' 2 . Al- 
though the mass scale at the maximum growth rate is com- 
paratively small, ~ ah 2 ~ 1 M©, larger mass scale per- 
turbation gradually grows in turn. As mentioned above, 
the gas density is higher in the outer part of the shell, 
then the unstable region extends from the outer edge to 
the inner part of the shell. By the age of Shl04, 0.7 Myr, 
the outer part of the shell dominated by CO molecule be- 
comes unstable. Therefore, the dense molecular fragments 
are formed around H II region and the triggered star for- 
mation may occur in the core of them, as DV03 shows. 

4.3. The Role of this Triggering Process 

Even if the molecular gas is successfully gathered around 
the H II region and the triggered star formation actually 
occurs, its impact on the global star formation is still un- 
certain. However, we can show its significance in a sim- 
ple argument. In our galaxy, most of stars are formed in 
clusters (e.g., Evans 1999), then the feedback effects from 
massive stars are ubiquitous for the global star formation. 
If every massive star gather the molecular gas of M s h on 
average and the triggered star formation occurs within the 
shell, the induced star formation rate SFRhii is estimated 
as 

SFRjjii ~ SFRo 
x I M ™ ^ 1 ( / (> 2QM ) \ / M sh \ / _e_\ 

yo.6M & J y 0.0006 J \iq a MqJ yo.ij 1 ' 

where SFR is the current total star formation rate in our 
galaxy, M av is the average stellar mass, / (> 20M Q ) is the 
number fraction of the massive star which creates relative 
large H II regions (> 20M©), and e is the star formation 
efficiency for the swept-up gas. The normalization values 
of M av - 0.6M© and / (> 20M Q ) ~ 0.0006 are calculated 
with the initial mass function (IMF) given by Miller & 
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Scalo (1979). Equation (2) shows that if the average mass 
of the gas shell is ~ 1O 4 M , this triggering process even 
alone can sustain the current galactic star formation rate. 
In our calculation, the molecular gas of ~ 10 4 M Q is accu- 
mulated around one massive star owing to the expansion 
of the H II region by the time of t <~ 1 Myr. There- 
fore, this mode of triggered star formation should be an 
important process. The efficiency of triggering should de- 
pend on the physical conditions, such as type of the central 
star, ambient density structure and so on. For example, 
many H II regions show the "blister-like" or "champagne 
flow" (Tenorio-Tagle 1979) features. In these cases, the 
IF rapidly erodes the parental cloud and only a part of 
the swept-up mass remains within the shell , which might 
result in smaller M s h (e.g., Whitworth 1979, Franco et al. 
1994). More realistic estimation including such cases will 
be explored in Paper II. 

5. CONCLUSION 

We have calculated the dynamical expansion of H II 
region and PDR around a massive star, focusing on the 
dense gas shell around H II region. Our results are simply 
summarized as follows, 



1 The dense gas shell is formed just in front of the IF, 
and this shell is initially atomic but gradually domi- 
nated with molecules. This is partly because the IF 
and SF overtake the preceding DFs and SF enters 
the molecular region, and partly because the refor- 
mation timescale within the shell becomes shorter 
than the dynamical time because of the high density 
of the shell. 

2 The calculation successfully reproduce the observed 
properties of Shl04. The mass of the shell and H II 
region, and the location of the IF at t ~ 0.7 Myr 
agree with the observation. The gravitational insta- 
bility occurs after the gas shell is dominated with 
the molecular gas. By the time of t ~ 0.7 Myr, the 
gravitationally unstable molecular shell is formed 
around H II region, as the observation shows. 

3 In our calculation, the molecular gas shell of the 
order of ~ 10 4 Mq is formed in 1 Myr around one 
massive star. With simple estimation, we show the 
triggering process owing to the expansion of H II 
region is important in context of Galactic star for- 
mation. 
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